
%Aspirin & 104 & 10933 & 11037 \\
%Placebo & 189 & 10845 & 11034 \\

aspi = [zeros(10933,1); ones(104,1)];
plac = [zeros(10845,1); ones(189,1)];
RR = (sum(aspi)/length(aspi))/(sum(plac)/length(plac));

BRR = [ ];
B=10000;
for b = 1:B
   baspi =  bootsample(aspi);
   bplac =  bootsample(plac);
BRR = [BRR (sum(baspi)/length(baspi))/(sum(bplac)/length(bplac))];
end

std(BRR)

figure(1)
hist(BRR, 80)
